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Abstract. We present a family of nonconforming vector finite elements of 
arbitrary order for problems posed on the space H(curl; Q) D H{div; Q), where 
n C M^. This result was first stated as a conjecture by Brenner and Sung 
in [l]. In contrast an extension of the same conjecture to domains of is 
disproved. 

Let Q he a domain of M'' where d £ {2,3}. As explained in [1] several problems 
involving the space (curl; r2)ni/ (div; $1), such as the cavity resonance problem and 
the acoustic fluid-structure interaction problem, can be solved using nonconforming 
finite element methods. In contrast conforming finite element methods cannot 
capture the solution of these problems under certain conditions. 

The accuracy of the approximate numerical solution of these problems can be 
improved if one uses finite elements which are not piecewise linear, but piecewise 
quadratic or of higher degree. For that purpose a quadratic nonconforming vector 
finite element for i?(curl; fJ) n i/(div;Sl) was introduced in [l], in the case of a 
two dimensional domain f2 C M^. The paper [I] also contains a conjecture which 
suggests a way of constructing nonconforming vector finite elements of arbitrary 
degree k for //(curl; ft) D H{div; ft), for domains of and of M.^. 

In order to state this conjecture and to formulate our results, we need to in- 
troduce some notations. We use boldfaced letters to represent vectors. The space 
of polynomials of total degree < fc in variables is denoted by Pfc(K''), and the 
space of homogeneous harmonic polynomials of degree fc in ci variables is denoted 
by 'Hk{^'^)- For each k >1 and d G {2, 3} we define a space Pfe^d of vector fields on 
M'' as follows 

(1) Vk,d [PkiM.'')]'' ® {VHk+2{M.'') © • • • ® VH2k{M.'')) 

For any triangle T if d = 2 (resp. tetrahedron T if d = 3) we consider a collection 
A4,d = ■^k d ^ -^k d of linear functional on Vk.d, which is defined as follows. The 
elements of ^{T) define the moments on T up to order fc — 2, for the d com- 
ponents of the vector fields. Denoting by Aq, ■ ■ • , the barycentric coordinates 
on the simplex T, and by ei, • ■ • , a fixed basis of R"^', these functionals have the 
expression 

M,,„(«) := / A;^''---A:^''(t;-eO, 
Jt 

for all 1 < i < 0? and all a = (ao, • ■ • , ad) G If^^ such -\- ■ ■ ■ + aa = k — 2. The 
elements of A/'|^(T) define the moments up to order k—low the d+1 edges (resp. 
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faces) of T, defined by Fj := {z G T : \j{z) = 0}, < j < d, again for the d 
components of the vector fields. These functionals have the expression 

for all 1 < i < d, aU < j < d, and all (3 = (J3o,--- ,/3d) G such that 

/3o + • ■ ■ + /?d — k — 1 and f3j = 0. We used in this definition the convention 
Xf = 0° = 1 on Fj. 

Brenner and Sung formulated in [T] a series of conjectures, which depend on the 
two parameters d G {2, 3} and fc > 1. 

, - Conj(/c, d) : For any T the elements of Vk,d are uniquely determined 
^ ' by the linear functionals in A4,d(r). 

The conjectures Conj(i,2) and Conj(l,3) are true and correspond to the noncon- 
forming Crouzeix-Raviart Pi vector finite element. It was established in [T] that 
Conj(2,2) is true, thus defining piecewise quadratic nonconforming vector finite 
elements in two space dimensions. 

The purpose of this paper is to establish the following result : 

Theorem. For any fc > 3 the conjecture Conj(fc, 2) holds. In contrast the conjecture 
Conj(2,3) does not hold. 

Our result therefore validates the construction of bi-dimensional vector finite 
elements of arbitrary degree proposed in [T] . On the contrary the three-dimensional 
quadratic vector finite element is invalid. Our result does not completely close the 
conjecture as the cases of three-dimensional vector finite elements of cubic or higher 
degree remain unsolved. 

It was established in [T] that for all d G {2, 3}, all fc > 1 and all T, one has 

dim Vk,d = i^^fk AT). 
Hence the conjecture Conj(fc,c?) is equivalent to the following property : 

, . For all T and all v G Vk,d, 

^ ' if {l{v) = for aU I G Mk.d{T)) then v = Q. 

In the first section of this paper we establish this property in the bi-dimensional 
case d = 2 and for an arbitrary fc > 1. In contrast we give in the second section a 
counter example in the three-dimensional case d = 3 and k — 2. 

1. Proof of the bi-dimensional result 

In this section the integer fc > 1 is arbitrary but fixed. If v — (vi, V2) G 'Pk,2 we 
remark that 

y-v = ^-^ePk.iiR^) and V..^|^ + |^GP._,(K^). 
0x1 0x2 0x1 0x2 

Our first lemma extends to degree fc an argument used in the initial paper [T]. 

Lemma. Let v G Vk.2- Let T he a triangle and let us assume that l{v) = for all 
l&MkAT)- Then 

(4) V-r> = Vxt> = 0. 
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Proof. We first notice tliat V-v and Vxv are polynomials of degree fc — 1, and 
that the components of V x (V x w) and of V(V • v) are polynomials of degree 
k — 2. In view of Green's theorem and the vanishing moments of v, we have 

{V X v){'V X v) dx — / {n X v){V X v) ds + / v ■ V x {V x v) dx — 
JdT Jt 

where n is the outer unit normal along dT. Similarly, we have 

/ (V • ■u)(V ■ v)dx ^ (n • ■u)(V ■ v) ds - v ■ V(V • v) dx = 0. 
Jt JdT Jt 

The results follow. O 

We now rephrase the conjecture ([2]) in terms of complex functions, and for that 
purpose we introduce some definitions. 

Definition. For any pair v = {vi, V2) of real valued functions we define a complex 
valued function P„ as follows 

P^{x + iy) vi{x,y) -iv2{x,y) for all {x,y) G R^. 

We now note that the equations \7 ■ v = Vx v ~ are equivalent to the Cauchy- 
Riemann equations of P^ , namely 

dx dy dy dx 

where 5R : C ^> M and S : C M respectively refer to the real and imaginary part. 
These equations characterize holomorphic functions. Let us introduce for all m > 1 
the space (D ,« of polynomials in the complex variable z = x + iy and of degree less 
or equal to m 

C™:-|p = f]^a,z'^; (ao, • • • , a,) G C"| . 

If V € 'Pk,2 satisfies l{v) = for all I G A4,2, then P^ satisfies the Cauchy-Riemann 
equations according to (|4]), and therefore P^ G C2fe-i. 

Definition. For any continuous function P : C ^ C and any zi, Z2 G C we define 

(5) /" P{zi+t{z2~ zi)){z2^ zi)dt^ f P{z)dz 
Jt=o Js 

where S C C is the oriented segment from zi to Z2 ■ 

Let S be an edge of a triangle T with endpoints {xi,yi) and {x2,y2) and let 
zi — xi+iyi and Z2 — X2+iy2 be their complex coordinates. Let v = {vi, V2) G Pk,2 
be such that l{v) — for all I G A/fc,2(T), and let Q{x + iy) :— Ri{x, y) + iR2{x, y) 
where Ri, R2 G Pfc_i(]R^) are arbitrary. Since v has vanishing moments up to order 
fc — 1 on the edges of T we have 

(6) Izi.Z2{PvQ) ^ {Z2- Zi) / {{viRl+V2R2)+i{viR2-V2Rl))\(x{t).y{t))d't 

Jt=o 

Where we used the notations x{t) := x\ + t(x2 — x\) and y{t) :— yi + t{y2 — yi). 
We now define a bilinear form which is related to our conjecture. 
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Definition. For all Z — (21,^2,23) G we define a bilinear form qz '■ C2/C-1 x 
(Cfe_i X Ck-i) C as follows 

qziP,{Qi,Q2)) -.^ iz.MPQi) + iz„zAPQ2)- 

Let T be a triangle and let zi — xi + iyi, Z2 = X2 + jj/2 and — + iy^ be the 
complex coordinates of the vertices of T. If 7; e 'Pk,2 is such that l{v) = for all 
I G J^k,2{T) then P„ G C2fc-i as previously noted. Furthermore, specializing ^ to 
polynomials Q G Cfc_i we obtain 

(7) qz{P.,iQi,Q2)) = for ah iQi,Q2) G Cfc_i x (Dk-i- 

The purpose of the rest of this section is to show that the bilinear form qz is 
nondegenerate. It then follows from Q that — and therefore that v — 
which concludes the proof of the conjecture Conj(fc, 2). 

We denote by B := {l,z, - ■ ■ , z^'^^^) the canonical basis of C 2k-i, and by B* := 
((l,0),(z,0)-- - ,{z''-\0), (0,1),-- - ,(0,z'=-i)) the canonical basis ofCfc_ixCfc_i. 
We denote by M(Z), or M(zi, Z2, z^), the matrix of qz in the basis B and B* . Hence 
for all 1 < i < 2fc and all 1 < j < fc we have 

Af(Z),,, - (z^-iz^-i) and M{Z),,j+k = I,,^,,{z'-^ z^-^) 

It follows that 

(8) Af (Z),, = - and M(Z),,,+, = .^"^^^ ■ 

The explicit expression of the matrix M{Z) = Af(zi, Z2, Z3), in the special case 
A; = 2, is given right after the end of the proof. 

Our next proposition gives an explicit expression of det M(Z), therefore showing 
that qz is non-degenerate. In the following Z always refers to the triplet of complex 
variables Z — (zi, Z2, Z3). 

Proposition. One has 

det M{Z) - a(zi - Z2)'=' (Z2 - Z3)'=' (Z3 - z,f 

where a = , V -'-"To, ?m > 0. Therej ^ore qz is non-degenrate whenever z\, Z2 and 
Z3 are pairwise distinct. 

Proof. We denote by S the collection of all permutations a of the set {1, • • • , 2fc}, 
and by s{<t) be the algebraic signature of such a permutation. We recall that 

2k 

(9) det M{Z) := ^ e{a) [] Af . 
For any permutation cr G § one has 

k k 

+ ^{3) - 1) + Y.^3 + ^(fc + 3) - 1) - 
3=1 i=i 

It follows from ([8]) that det M{Z) is a homogeneous polynomial in the variables 
^ij Z2,zz and of degree 3fc^. We also note for future use that 

k 

(10) E(j' + ^ovi)>fc' 
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with equality if and only if a leaves invariant the sets {1, • • • , fc} and {k+1, • • • , 2k}. 
For any c G (D we define two 2k x 2k triangular matrices P{c) and P*{c) associated 
with the following changes of basis on C2fc-i and (Ck-i x (C k-i respectively 

P{c)B = ( 1, z + c,-- - , (z + c)2'=-i) 
P*{c)B* - ( (1,0),.-. , ((z + cf-\0), (0,1),--- , (0,(z + cf-i)) 

One easily sees that the matrices P{c) and P*{c) are upper-triangular and have 
ones on the diagonal, hence det P(c) = detP*(c) = 1. 
Since 

iz,+c,z^+c{z'z^) = h,,z^{{z + cy{z + cy) 

we obtain 

M{zi + c, Z2+ c, Z3 + c) = P{cfM{Z)P*{c). 
Recalling that det P{c) = det P*{c) = 1, and choosing c — —zi, we obtain 
detM(0, Z2-ZI, zg - zi) = det Af(Z). 

The explicit expression of the matrix M(0, Z2 — zi, Z3 — zi), in the special case 
= 2, is given right after the end of the proof. It follows from (jS]), ^ and (fTO|) 
that the polynomial det M{Z) is a multiple of (z2 — zi)'' . Similarly, det M{Z) is a 
multiple of (za — zi)*^ . 

Substracting column fc + i from column i, for all 1 < i < fc, we find that 

detM(z3,Z2,zi) = {-i f det M{Z) 

and therefore detAf(Z) is also a multiple of (Z3 — Z2)''' . Since detM(Z) is a 
polynomial of degree 3fc^ in the complex variables zi, Z2, Z3, and since (zi — Z2)'^ , 
(z2 — z^)'^ and (Z3 — zi)*^ have no common factors, there exists a constant a £ (D 
such that 

det MiZ) = a (zi - Z2f (za - Z3)'=' (zg - zif . 
In order to compute the constant a, and to show that a ^ 0, we remark that it 
is the coefficient of z'^ in the polynomial detM(0,z,l) = a(-z)'= (z - 1)'= . The 
explicit expression of this matrix, in the special case fc = 2, is given right after the 
end of the proof. 

The contribution of a permutation ct G § to det M(0, z, 1) is a monomial which 
has degree k'^ if and only if ([TU]) is an equality. Denoting by S* the collection of 
permutations of the set {!,••• , k}, we obtain that detM(0, z, 1) equals 

Hence using © 

detAf(0,z,l) = z'='det ( ^ ) det ( ) + ©(z'^'+i) 

V « + i - 1 / i<tj<k \i + j + k-lj i<ij<k 

This expression gives the value of a as the product of two Cauchy determinants, 
which can be computed using the formula, established in [5] §1.1.3, 



det 



n ("i - aj) n {bi - bj 

1 \ l<i<j<k l<i<j<k 



a« + ^j7i<,.j<fe n (ai + bj) 

~ ' ~ l<i.j<k 
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This concludes the computation of det M{Z). o 

The following explicit expressions of the matrix M, if = 2 and z,Zi,Z2,Z3 e C , 
may be useful to the reader 

/ 2 2 2 2 \ 

/ Z2-Z1 ^2-^ - ^2_lL \ 



M{zi,Z2,Z3) = 



M{0, Z2 - Zi, Z3 - Zi) 



M{Q,z,l) 



Z3 - Zi 

2 2 



/ -22 - -21 

(Z2-Zl)^ 

2 

(z2-z,y' 



5 

(^2-^l)^ 

2 

(Z2-Zl)^ 

(Z2-Zl)^ 

4 ^ 

(^2-^l)" 



zz - Zi 

(Z3-Zlf 

2 ,, 

{z3-zir 

4 



{Z3-Zif 

2 

(Z3-.^1)" 



-2 .3 , , 



2 

2 3 



^4 5 



1 1 

2 3 
1 1 

3 4 
1 1 

4 5 



2. A COUNTER EXAMPLE IN THREE SPACE DIMENSIONS 

Let To be the simplex of vertices (0, 0, 0), (1, 0, 0), (0, 1, 0), (0, 0, 1) and let Pq 
be the harmonic polynomial of degree 4 



Po 



3x 



lOx^ - 15.T'* + 3y - 18xy - Ihx^y + SOx^y - l^xy"^ + 45a:^y2 



+2,Qxy^ - Ihy"^ + 3z - l%xz - Ibx^z + SOx^z - 18yz + 2A{)xyz - ld,Qx^yz 
-15y^z - ISOxy^z + 30?/z - 15x2^ + A^x^z^ - 15yz^ - ISOxyz^ + Aby'^z^ 
+10z^ + 30xz^ + 30yz''^ - 15z*. 
We define 

uo := VPo e ^2,3- 

One can easily check using a formal computing program that all the linear func- 
tionals in A/2,3(To) vanish on uq, which shows that the conjecture Conj(2,3), on 
quadratic vector fields in three dimensions, is not valid. The interested reader can 
download on the website arxiv.org, jointly to the preprint of this paper, a file that 
contains these verifications. 

Finding a quadratic vector finite element for iJ(curl; fi) ni?(div; 0.) in dimension 
3 thus remains an open question. Let us finally mention that, up to a multiplicative 
constant, uq is the only element of 7^2,3 on which all the linear functionals A/'2,3(To) 
vanish. 
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